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I.  INTRODUCTION 


The  object  of  this  study  is  to  develop  an  efficient  numerical  solu¬ 
tion  for  the  high-frequency  near  zone  radiation  patterns  of  a  spheriod- 
mounted  antenna.  The  Geometrical  Theory  of  Diffraction  (GTD)  is  em¬ 
ployed  for  the  pattern  calculation  of  which  the  creeping  wave  [1]  solu¬ 
tion  in  the  shadow  region  is  of  particular  interest  here.  According 
to  the  generalized  Fermat's  principle,  a  ray  emanating  from  a  source, 
which  is  located  on  the  surface,  follows  a  geodesic  path  on  the  surface 
and  continually  sheds  energy  into  the  shadow  region.  Such  a  creeping 
wave  mechanism  is  illustrated  in  Figure  1(a),  from  which  it  can  be  seen 
that  a  ray  traverses  from  the  source  point  Q'  to  the  diffraction  point 
Q,  and  then  propagates  along  the  geodesic  tangent  at  Q  toward  the  ob¬ 
servation  point  P^.  Since  most  parameters  involved  in  the  pattern  cal¬ 
culation  are  related  to  the  ray  path,  the  geodesic  solution  on  the  surface 
for  a  near  zone  field  point  constitutes  a  major  portion  of  this  study. 


(a)  FIELD  POINT  IN  b  *  b^' 

SHADOW  REGION  FIELD  POINT  IN 


LIT  REGION 

Figure  1.  Ray  paths  in  the  shadow  and  lit  regions. 
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The  spheroid  is  studied  because  of  its  three  dimensional  nature 
in  that  it  can  be  used  to  simulate  an  aircraft  or  missile  fuselage. 

For  a  flush-mounted  airborne  antenna,  the  fuselage  has  a  dominant  effect 
on  the  resulting  radiation  pattern.  Previously  [zl  the  fuselage  was 
simulated  by  an  elliptic  cylinder;  however,  the  cylinder  model  could 
not  predict  the  pattern  close  to  the  nose  or  tail  sectors,  where  the 
deviation  from  the  physical  situation  becomes  very  prominent.  Previous 
attempts  [2 1,  using  the  spheroid  model,  employed  calculus  of  variations 
or  tensor  analysis  to  calculate  the  geodesic  paths.  The  resulting  al¬ 
gorithms  were  very  time-consuming  and  virtually  impractical  to  use. 
Therefore,  a  perturbation  method  is  presented  here  in  an  attempt  to 
combine  the  simplicity  of  the  cylinder  or  cone  solution  with  the  gen- 
e>^ality  of  the  prolate  spheroid. 

Before  examining  our  perturbation  approach,  it  is  appropriate  to 
study  the  basic  radiation  mechanisms  in  some  detail  in  that  they  suggest 
the  present  solution.  The  ray  analysis  in  the  lit  region  can  be  rather 
straightforwardly  analyzed  for  any  convex  surface;  however,  the  geodesics 
which  are  necessary  to  treat  the  shadow  region  dictate  what  surfaces 
can  be  efficiently  analyzed.  In  the  shadow  region,  the  energy  propagates 
outward  from  the  source  along  the  geodesic  paths.  As  the  energy  flows 
around  the  surface,  it  is  continuously  diffracted  along  the  geodesic 
tangent  toward  the  field  point  such  that  the  significant  effect  of  the 
surface  is  associated  with  a  region  around  the  source.  In  fact,  for 
a  prolate  spheroid,  the  significant  portion  of  the  surface,  which  is 
associated  with  the  dominant  energy,  may  look  as  shown  in  Figure  2. 

Note  that  this  region  can  be  specified  by  following  the  various  geodesic 
paths  until  the  radiation  level  along  a  given  path  becomes  insignificant, 
i.e.,  more  than  40  dB  below  the  pattern  maximum.  With  this  in  mind, 
it  is  clear  that  one  could  represent  the  prolate  spheroid  by  a  structure 
which  simulates  the  circular  cross-section  completely;  however,  the 
profile  (i.e.,  the  elliptic  curve)  could  be  approximated  by  a  simpler 
shape  in  that  the  significant  energy  region  does  not  cover  a  large  portion 
of  the  profile  shape.  Consequently,  a  right  circular  cone  perturbation 
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is  chosen  here  to  simulate  the  spheroid  which  in  turn  can  be  used  to 
model  a  fuselage.  This  perturbation  model  is  illustrated  in  Figure 
3  for  a  source  located  near  one  end  of  the  spheroid.  Note  that  if  the 
antenna  is  placed  at  the  center  of  the  spheroid,  the  cone  actually  be¬ 
comes  a  right  circular  cylinder. 

In  that  the  cone  is  a  developed  surface,  one  can  unfold  the  cone 
such  that  a  planar  structure  results.  The  geodesics  associated  with 
the  cone  are,  then,  straight  lines  on  this  planar  structure.  In  order 
to  allow  for  a  geodesic  solution  between  the  simplicity  of  the  cone 
and  the  rigor  of  the  spheroid,  one  can  perturb  the  cone  by  bending  it 
along  its  generator  as  illustrated  in  Figure  3(b).  In  that  a  perturba¬ 
tion  technique  is  employed,  the  geodesic  paths  for  the  cone  are  simply 
modified  such  that  the  solution  for  the  spheroid  is  basically  straight¬ 
forward  and  is  the  basis  for  this  report.  It  is  obvious  that  one 
cannot  use  this  perturbation  if  significant  energy  propagates  great 
distances  along  the  surface.  However,  as  discussed  above,  the  energy 
which  propagates  great  distances  along  the  spheroid  generator  becomes 
insignificant  in  magnitude  such  that  one  need  not  solve  for  the  true 
geodesic  paths  outside  the  region  shown  in  Figure  3.  The  simplicity 
of  these  perturbed  geodesic  paths  allows  one  to  very  efficiently  determine 
the  significant  ray  paths  on  the  prolate  spheroid. 

It  is  noted  that,  since  GTD  is  a  high-frequency  method,  the  lower 
frequency  limit  of  this  solution  is  dictated  by  the  electrical  dimen¬ 
sions  of  the  spheroid,  i.e.,  the  semi-minor  axis  of  the  spheroid  is 
required  to  be  at  least  a  wavelength. 

The  curved  surface  diffraction  solutions  employed  here  are  briefly 
discussed  in  Section  II.  The  parameters  needed  for  this  GTD  solution 
are  described  in  Section  III  in  terms  of  the  cylinder  and  cone  pertur¬ 
bations  for  near  field  used  to  simulate  the  geodesic  paths  on  the  spheroid. 
Computational  results  are  presented  in  Section  IV.  Finally,  a  summary 
of  this  study  is  made  in  Section  V. 
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II.  THEORETICAL  BACKGROUND 


A.  Introduction 

The  radiated  field  of  a  spheroid-mounted  antenna  is  analyzed  using 
the  Geometrical  Theory  of  Diffraction  (GTD).  The  surface  is  assumed 
to  be  perfectly  conducting,  and  the  surrounding  medium  is  free  space. 

An  exp  (jtot)  time  dependence  is  understood  and  suppressed  in  the  fol¬ 
lowing  formulations. 

Consider  an  infinitesimal,  magnetic  current  moment  or  an 

electric  current  moment  dPg(Q‘)  located  on  a  perfectly  conducting  con¬ 
vex  surface  as  shown  in  Figure  1;  the  sources 


flPm(Q')  =  E(q')xn'da'  and  dPg(Q')  =  I(ll')difn' 


pertain  to  the  aperture  and  monopole  type  excitations  with 
r{Q')  =  electric  field  at  Q', 
n'  =  outward  unit  surface  normal  at  Q', 
da*  =  area  element  at  Q', 

1(4')  =  electric  current  distribution  on  the  monopole,  and 
£'  =  distance  parameter  along  the  monopole. 
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According  to  geometrical  optics,  the  space  surrounding  the  source  is 
divided  into  an  illuminated  and  shadow  region  by  a  plane  tangent  to 
the  surface  at  Q'.  This  plane  is  referred  to  as  a  shadow  boundary.  Since 
the  field  in  the  deep  lit  region  is  essentially  that  obtained  from  ge¬ 
ometrical  optics,  and  the  field  in  the  deep  shadow  region  is  relatively 
weak,  the  solution  for  the  transition  region  adjacent  to  the  shadow 
boundary  is  of  more  interest  and  discussed  below. 


Shadow  Region 


The  creeping  wave  mechanism  in  the  shadow  region  is  illustrated 
in  Figure  4.  From  the  generalized  Fermat's  principle,  a  ray  emanating 
from  the  source  dPj^(Q')  at  Q'  traverses  a  geodesic  path  Q'Q  on  the  sur¬ 
face,  and  propagates  along  the  geodesic  tangent  at  Q  toward  the  field 
point  P^.  The  field  dEjg  at  can  be  expressed  in  terms  of  the  field 

at  a  reference  point  P^  by  [1] 

-jks 

<l%(Po)  7T— J—  e  °  (1) 

J  (pi+Sq){P2+Sq) 

where 

P  ^  and  are  the  principal  radii  of  curvatures  of  the  wavetVc.it 

p  Q 

and  oDn”  ,m“'^,...3  are  the  higher  order  terms. 


From  Figure  4,  it  is  seen  that  if  the  reference  point  P  is  moved  to 

H  H  ^ 

the  curved  surface  diffraction  point  Q,  then  P^-^,  P^  and  s^-^s. 
Since  dE^(P^)  is  independent  of  the  reference  point  P^,  it  follows  that 


Pq-qJ^  dEg(PQ)  =  l?(Q\Q)  ; 

pJ-K) 


(2) 
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(8) 


and 


(9) 


9(0=7^/" 


i  exp(-j2iT/3) 


9(0  -  ^! 


exp(-jT£ 


yir  "  exp(-j2Tr/3) 


which  are  known  as  the  acoustic  hard  and  soft  Fock  functions.  The  Fock 
type  Airy  function  is  given  by 


Wo(t)=—  r  dt-exp(Tt-t^/3)  (11) 

00  exp(j2u/3) 

and  w^'(  t)  is  the  derivative  of  W2(t)  with  respect  to  t.  The  Fock  para¬ 
meter  ^  for  the  shadow  region  is  given  by  Cl] 


_  _  jn  m(t’) 

f,  -  J  at' 


kP„(f ) 


ni(t')  =  — ^ 


Here  Pg(t')  is  the  surface  radius  of  curvature  along  the  ray  path  at 
t'.  The  width  of  the  surface  ray  tube  at  Q,  dn(Q),  is  given  by 


dn(Q)  =  p^di|)  . 


The  parameters  and  t  are  defined  as  the  free  space  wave  impedance 
and  geodesic  arc  length  from  Q'  to  Q,  respectively. 


Combining  Equations  (3)-(14),  the  n  and  b  directed  components  of 
dr(^i(l’^)  are  given  by  [l  1 


dPm(Q')  case: 


dE>s) 


s(p  +s) 


e-jks+o[m‘2] 


[(dPm.b)T^S+(dPm.t')SIe-J''^^a^] 


0  /  1 


s(Pc+s) 


[m'^,m'^] 


(b)  dPg(Q')  case: 


dE"(Ps) 


-jkZ  p  (Q')  /— 

dPg(Q-)He  1  J  dWs(p^-.s) 


e''^'^^+0[m'^] 


=  TT^ 


;(Q  )TQSe  ^pjQ)  J  J  dt  Js(p^+s) 


e'J''^+0[m’^] 


where  To=T(Q' )pg(Q' )  with  T(Q')  being  the  surface  torsic,.  at  the  source 
location  (refer  to  Table  1). 
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^  T 


C.  Lit  Region 


From  geometrical  optics,  the  source  dP|g(Q')  at  Q'  excites  waves 
which  propagate  along  straight  line  ray  paths  from  the  source  to  field 
point  in  the  lit  region.  As  shown  in  Figure  5,  the  field  cirg{P|^  at 
point  Pj^  is  expressed  by 


dF@(PL)  -  dE@(f,)  ^  -  e"'  '°-o[m-2.m-3]  (19) 

J  (Pi%)(p2+So)  L  J 

Since  Q'  is  the  only  caustic  of  the  incident  rays,  the  principal  radii 

1  4  v/' 

of  curvature  Pj  and  associated  with  the  incident  wavefront  at  P^  are 

identical,  i.e.,  pj  =  p’.  Furthermore,  dEig(P|^)  is  independent  of 

the  reference  point  If  P  is  chosen  to  be  at  O',  it  follows  that 

0  0 


i  ~  ^ 

^im  p  dEg  (Po)  =  i-g 


should  exist.  Thus,  Lg  can  be  related  to  dP|g(Q')  by  [1] 

I 

=  dPg(Q')  • 


Equations  (19)-(21)  are,  then,  combined  to  yield 
dEig(PL)  ^  dPig(Q')  •  Tg-  — - —  +  0[mj^  ,mj^f...  ] 


■  if  Il>inAnibB*b'bjC*tj.,D] 


— aO  ^  ^  /V 

In*nM+n*bN] 


with  A,  B,  C,  D,  M,  and  N  are  defined  in  Table  2.  Note  that  dt^(P|  ) 

/N  e  L 

is  decoupled  into  n  and  b  components  as  follows: 


13 


(a)  case: 


=  ;^[(dP^*b*)(H*’+T^rcosoM  +  {cr^•^  )T/cose’ 1 


r  -2i 

V-"  oK  ^ 


"EmCl)  =  if [  («'„-b')V(dP„.t‘’)(S«’-T2Fcos26')] 


^“jks  A  o 

7—  *  OCm;  ,m;  ] 


(b)  ^g(Q')  case: 


-jkZ 


"  “4^  clPg(Q')sine’ IVl^+T^fcose’ ]^-g — ^  0  3 


-jkZ, 


-jks 


clP3(Q')sine\FV-0  6’-/3 


where 


H‘=  9(5()e 


-j4'3 
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and 


=  -nij^(Q')cose’ 


(31) 


m  (O')  =  - "!(.Q‘.) _ 

(UTW0i)l/3 


(32) 


The  angle  0^  is  defined  by  fi'.§  =  cos0’  as  shown  in  Figure  5,  Also, 


F  =  S^-H^ose^ 
1+T^cos^  e’ 


(33) 


as  defined  in  Table  IT. 


D,  Pattern  Factors 

The  solutions  for  short  magnetic  or  electric  dipoles  have  been 
given  in  parts  B  and  C.  One  approach  to  analyze  an  extended  aperture 
or  linear  antenna  problem  is  to  integrate  the  above  solutions  over  the 
source  distribution,  if  it  is  known.  This  is  an  application  of  the 
superposition  theorem,  and  one  approximates  the  source  distribution 
by  an  array  of  short  magnetic  (or  electric)  dipoles  on  the  conducting 
surface.  This  is  an  accurate  solution,  however,  rather  tedious.  A 
more  efficient  approach  is  to  modify  dP^(Q')  such  that 


(a)  in  the  shadow  region: 


cos 

(r  c.-f'))' 

sin 

' 

1 - 

• 

^e  ~  n'[l-cos(kL)j 


(35) 
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(b)  in  the  lit  region: 


cos 

' 

sin  ^  sine^P^.b') 

1'. 

•p  sineVP^*b')  . 

(36) 


F 

e 


cos(kLn'«s)-cos(kL) 

l-(n'«s)^ 


(37) 


Here  =  unit  vector  in  the  direction  of  magnetic  current  moment, 

A,B  =  the  length  of  the  short  and  long  sides  of  the  slot,  and 
L  =  the  length  of  the  monopole. 

It  is  noted  that  L  is  not  to  exceed  a  quarter  wavelength  for  the  solu¬ 
tions  to  be  valid. 

III.  NUMERICAL  TECHNIQUE  AND  PERTURBATION  METHOD 
A.  Introduction 


It  is  seen  that,  for  an  antenna  mounted  on  a  spheroid,  the  solu¬ 
tion  in  the  lit  region  is  relatively  straightforward.  On  the  other 
hand,  the  geodesic  paths  associated  with  the  GTD  solution  in  the  shadow 
region  are  extremely  complex.  Previous  studies  [2]  employed  calculus 
of  variations  or  tensor  analysis  to  calculate  the  geodesic  paths  which 
resulted  in  very  inefficient  solutions.  An  efficient  numerical  approach 
is  examined  in  this  chapter  with  the  spheroid  simulated  by  a  perturbed 
cone  or  cylinder  model.  In  that  the  cone  and  cylinder  are  developed  sur¬ 
faces  for  which  the  geodesics  are  efficient  to  solve.  Given  a  radiation 
direction  (6^.,(t)^),  one  can  find  the  final  diffraction  point  (eQ,ij)^) 
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by  following  the  geodesic  path,  step  by  step,  until  the  geodesic  tangent 
coincides  with  This  is  a  rather  tedious  and  time  consuming 

process  if  applied  for  each  new  radiation  direction.  Considering  t*>at  a 
new  radiation  direction  does  not  deviate  greatly  from  the  previous  di¬ 
rection,  one  should  be  able  to  develop  a  solution  which  uses  the  prop¬ 
erties  of  the  surface  and  the  previous  geodesic  path  to  find  the  new 
diffraction  point.  Such  an  approach  is  applied  here  to  make  this  solu¬ 
tion  as  efficient  as  possible. 

Since  the  field  decays  exponentially  along  the  ray  path  on  the  sur¬ 
face,  it  is  assumed  that  only  one  or  possibly  two  dominant  rays  exist 
in  the  problems  treated. 

B.  Spheroid  Surface  Parameters 

Using  the  spheroid  geometry  shown  in  Figure  6,  the  surface  is  de¬ 
fined  by 


1l(0,(j))  =  R(0)sin0cos4»x+R(0)sin0sin(t)y  +  R(0)cos02 


(38) 


with 


R(0) 


_ab_ 


J 


9  9  7  7 

a  COS  0+b  sin  e 


ro) 


or 


R(v,(())  =  a  cos\coS(|)x  +  a  cosvsin(t)y  +  b  sinvz 


(40) 


with 


tanv  .^-1® 
a  SI n6 


(41) 
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The  "v"  parameter  is  introduced  because  of  its  convenience  in  analyzing 
elliptic  geometries. 

Considering  a  ray  which  propagates  along  a  geodesic  path  Q'Q  on 
the  spheroidal  surface  as  shown  in  Figure  7  the  three  unit  vectors  t, 

/N  ^ 

n  and  b  are,  as  defined  earlier,  the  geodesic  tangent,  outward  surface 
normal  and  binormal  at  any  point  along  the  geodesic  path.  The  outward 
surface  unit  normal  (n)  is  obtained  from 


n  = 

lR(|)xRvl 


where 


(42) 


V  =  -a  sin\)COS(t)X  -  a  sinvsin(j>y  +  b  cosvz 


and 


3R  _ 


94> 


=  R(|)  =  -a  cosvsin(|)x  +  a  cosvcos<|)y 


Then,  the  unit  normal  vector  is  given  by 


n  =  b  cosvcos(t)X  +  b  cosvsin<}>y  +  a  sinvz 

p  2  2~ 

J  a  sin  V  +  b  cos  v 


(43) 


The  normal  curvatures  on  the  surface  are  evaluated  by  introducing 
the  first  and  second  fundamental  forms  of  differential  geometry  [3I 
such  that 


-  Ld4>^  +  2Md<l)dv  +  Ndy^ 
Ed4i^  +  2Fd(t)dv  +  Gdv^ 


(44) 
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a> 


Figure  7.  Geodesic  path  on  a  spheroid. 
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where 


1-  -  RA<^•N.  M  =  R^^.N.  and  N  = 

^  ^  ^  =  W' 

and 

> 

N  =  -n  . 

It  can  be  shown  that 

=  -a  cos\)Cos<bx  -  a  cos  sinchy  , 

<p<p 

=  a  sinvSln(j)X  -  a  sinvcos<j)y  ,  and 
R^^  =  -a  cosvcoS(|)X  -  a  cosvsirvjjy  -  b  sin^z. 
After  some  algebraic  manipulation,  one  obtains 
L  =  ab  cos^v/y?sin^v  +  b^cos^v  . 

M  =  0  . 

N  =  ab/Ja^sin^v  +  b^cos^v  . 

?  7 

E  =  a^cos  V  > 

F  =  0 .  and 

7  7  7  7 

G  =  a  Sin  V  +  b  cos  v  • 
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The  two  principal  curvatures  on  the  surface  are,  then,  given  by 


=  L/E  = 


19  2  9  2 

ava  sin  V  +  b  cos  v 


Ko  ~  N/6  =  . — 


’'(a'^sin  V  +  b^cos  v)' 


which  can,  also,  be  expressed  as 


_  c  /a^cos^e  +  b^sin^e 
^  ^7  a^cos^e  +  b^sin^e 


-  ab(a^cos^e  +  b^sin^e)^^^  /.^v 

It  is  noticed  that  Rj  =  ^  and  R2  =  ^  as  found  in  Table  1.  Since  F*M=0, 
the  directions  and  ^  ire  in  the  principal  surface  directions  [3] 

4*  y  .s  r  r  L  j 

which  are  expressed  by  Tj,  and  are  associated  with  kj  and  <2* 
spectively,  as  shown  in  Figure  7.  If  3  denotes  the  angle  between  t 

A 

and  Tj,  then 

A  A 

t  =  Tjcoss  +  T2Sin3  .  (49) 

From  Euler's  theorem,  the  normal  curvature  along  the  geodesic  path  is 
specified  by 

2  2 

Kg  '  Kjcos  B  +  K2Sin  6  (50) 


with  the  radius  of  curvature  (pg)  being  1/Kg. 

The  torsion  term  (T^)  introduced  in  Section  11(B)  is  given  by 
^  Pg 

where  the  surface  torsion  is  given  by 

T=^^(^-Ki)  (52) 

with  Kj  and  being  defined  in  Equaitons  (45)-(48). 

C.  Surface  Geodesics 


The  geodesics  on  the  cylinder,  cone  and  spheroid  are  examined  in 
this  section. 

(a)  Cylinder  case: 

The  circular  cylinder  geometry  used  for  this  study  is  shown  in 
Figure  8(a).  Since  the  cylinder  is  a  developed  surface,  the  geodesic 
path  Q'Q  is  a  straight  line  on  the  unfolded  planar  surface.  As  shown 
in  Figure  8(b),  the  geodesic  unit  tangent  t  is  given  by 


t  =  (j)COSY  +  zsiny  (53) 


where 


<|)  =  -sin(|)x+cos4)y. 


(54) 


It  is  noticed  that  41  and  z  are  the  two  principal  directions  on  the  cyl¬ 
inder  surface.  For  a  given  geodesic  Q'Q,  one  can  see  that  y  is  a  con¬ 
stant  along  the  geodesic  path. 
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(b)  UNFOLDED  PLANAR  SURFACE 


Figure  8.  Geodesic  path  on  a  developed  cylinder. 


(b)  Cone  case: 


Consider  a  ray  propagates  along  a  geodesic  path  Q'Q  on  the  cone 
surface  as  shown  in  Figure  9(a)  which  is  a  straight  line  on  the  unfolded 
planar  surface  as  shown  in  Figure  9(b)  since  the  cone  is  a  developed 
surface.  Applying  the  same  technique  as  in  Section  B,  t^  and  ^  are 
the  two  principal  directions  on  the  surface.  The  geodesic  unit  tangent 
t  is,  then,  represented  by 

/V 

t  =  (|cos8+tgSinB  (55) 

where 

/N  Ai  /\ 


=  -XgSin6+2cos6  , 

(56) 

=  cos<|>x+sin<t'y. 

(57) 

and  6  is  the  half  cone  angle  as  shown  in  Figure  9(a).  Note  that  3  is 
no  longer  a  constant  along  the  geodesic  path  Q'Q. 

(c)  Spheroid  case: 

Recall  that  the  prolate  spheroid  surface  is  defined  by 

A  A 

R(v,<j))  =  a  cosvcos(tix  +  a  cosvsin(|)y  +  b  sinvz 


with 


tanv  = 


b  cosQ 
a  sin<|>  * 
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— — 


From  Section  B,  it  was  found  that  the  principal  directions  were  given 
by  and  As  seen  in  Figure  10(a)  and  (b)  the  unit  vectors  in  the 
principal  directions  can  be  represented  by 

tr  =  -xsin4i  +  ycoS(J)  (58) 

t^  =  R  /l^  I  =  -x„sin6+zcos6 
e  v  *  v'  e 

A 

-a  sinvx^  +  b  cosvz 

=  - - - ^  (59) 

Ja^sin^v+b^cos^v 

wi  th 

Xg  =  cos4»x+sin4iy  . 

A  A 

It  is  noticed  that  t  is  merely  <j).  If  B  denotes  the 
and  the  geodesic  unit  tangent  t,  then 

A  A  A. 

t  =  tj.cosB+tgSin6  ,  (61) 

which  is  identical  to  the  form  used  for  the  cylinder  and  cone  geodesics. 
This  suggests  that  one  might  be  able  to  develop  a  perturbation  solution 
which  gives  a  simplified  form  for  B  on  a  spheroid  using  the  cylinder 
or  cone  expressions  for  g. 

D.  Cylinder  Perturbation  for  Near  Field 

In  order  to  solve  for  the  geodesic  paths  on  a  prolate  spheroid 
for  the  near  zone,  the  cylinder  perturbation  technique  is  used  when  the 
source  is  located  at  or  near  0^^9O°.  As  shown  in  Figure  11,  the  cylinder 
is  perturbed  by  bending  it  around  the  spheroid,  and  the  geodesic  paths 
are  to  be  resolved  in  this  perturbed  cylinder  model. 


(60) 

/s 

angle  between  t^ 
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Figure  10(b).  Elliptic  profile. 

Recalling  that  y  is  a  constant  along  a  given  geodesic  path  Q'Q 

on  the  perturbed  cylinder  in  Figure  11,  one  obtains  a  geodesic  equation 

given  by 

Se 

tanY=5-  (62) 

r 

with  =  3(1^,  where  <t^  is  measured  from  the  source  to  the  diffraction 
point.  The  elliptic  cross-section  in  the  x^-z  plane  is  described  by 

Xg  =  a  cosVg  and 

z  =  b  sinM 

e 

therefore,  the  arc  length  is  given  by 

'’e  r — - - 

Sg  =  /  va^sin^v  +  b^cos^v  dv 
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Thus,  Equation  (62)  becomes 


tan  Y  = 


/  Ja^sin^ V  +  b^cos^ v  d v 


>4^. 


(63) 


with  -Tr<Y<Tf>  “  n/2£(ty.,  Vg£Ti/2. 


Assuming  the  diffraction  point  Q  (a  cosVgCOs4)j.,  a  cosVgSincty,  b  sinvg) 
and  the  receiver  point  (R^sin0^coS()(^,  R^sine^sin(|i[^.,  R|.cose^),  then 
at  the  diffraction  point  Q,  the  radiation  direction  (©^.(Ixj.)  should  co¬ 
incide  with  the  geodesic  tangent  t  given  in  Equation  (61).  Thus, 


t  =  X  t^  +  y  t^  +  z  t^ 

y\  ^ 

=  t^  COSY  +  tg  siny 


(64) 


where. 


R^sin0^cos<()^  -  a  cosVg  cos(|)j. 


R^sin0^.sin<t)^  -  a  cosVg  sin<j(^ 


R.  cos0^  -  b  sinv 
t^  =  — - 6 - ^  .  and 


=  R^-2R^[a  sin0^cosvgcos(<|)^-<j)^)+b  cosO^sinVg^ 


+  a^cos^  +  b^sin^ 
e  i 
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tquating  the  x-,  y-,  and  z-components  respectively,  one  obtains 


-Sintf^COSy 


a  s1nvgCOS{jysin  y  R^sin0^coS(j)^  -  a  cos^gCoscj)^ 


sin  V  +  b  cos  v 

C  €? 


(65) 


COS(j)  cosy 


a  sinVgSin(|)^sinY  R^sin9^sin(j)^  -  a  cosVgSincj)^ 

f~P  ^  ^  p  D 

Ja  sin  +  b  cos  V 
e  e 


(66) 


b  cosVgSin  y 


R^  cosG^  -  b  siny;^ 


J  a  ^sin^  +b^cos^v. 


(67) 


From  equations  (65)  and  (66),  one  obtains 


cos  Y 


R^.sin9^sin(<|»^-<))^) 


(68) 


and 


a  sinVgSiny  R^sin0^cos(())^-())^)-a  cosvg 

I  2  .  ?  ?  ^ 

Ja  sin  V  +b  cos^Vq 


(69) 


Equation  (69)  is  divided  by  (67)  to  give 


-a  sinvg  R^sino^cos(  (|)n.-()y) 
i)  cosVg  ”  R^cose^-b  sfn^ 


(70) 
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From  Equations  (62),  (67)  and  (68),  one  obtains 


b  cosVgSg  R^cose^-b  sinvg 

loo  99  "  RtS^'’6tSin((j)t-(l>^) 

a4)Ja^sin^Vg+b^cos^Vg 


From  Equations  (70)  and  (71),  one  obtains 


(71) 


f(R^,0^,<{)^,Vg,<|)^)  =  b  sine^cos((|)^-<|)^)cosvg+a  cose^sinvg 
ab 

=0  ,  and  (72) 


g(R^,0t,4)t,Vg,(t)^)  =  b  $gSin0^sin((|)^-^,^)cosvg-a4,^Ja^sin^Vg+b^cos^Vg 


COS0^  + 


ab(t) 


rv 


?  2  2  2 
a  sin  V  +b  cos  v^sinv^ 
e  e  e 


R* 


=  0 


(73) 


Provided  that  one  has  obtained  a  diffraction  point  (Vg.ct)^)  for  a  receiver 
location  (R^,0^,<i)^) ,  a  numerical  technique  can  now  be  developed  from 
Equations  (72)  and  (73)  to  solve  for  (Vg+AVg,  '|)^+A())^)  associated  with 
a  new  receiver  location  (R^+AR^,  Q^+A0^(>^+A(|)^) .  Assuming  that  the 
ith  set  of  (9t’'^’t’'’e’***e^  first  known  to  satisfy  f^.=g^=0,  or  at  least 
approximately  so,  the  next  set  (R|.+AR^,0|.+A9|.,((>^+A(|)^,Vg+AVg,(l)^+A<j)^) 
is  obtained  by  enforcing  such  that 


and 


9i+l 


=  Qi+gR^AR^  +  g^^AGt  + 


"  9,  AVe 


+  9 
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In  matrix  form,  it  is  given  by 


I 

I 


I 


Av^ 

Ve  <D^ 

UVg 

A4)r  _ 

L 

-fi  -  aR^  -  ■  ^(t)^^‘'’t 


L-Qi  ARt-9e^^®t  -\^h 


fy  =  -b  sinS^cos (4)^-4)^) sinVg  +  a  cosO^cosVg 
e 

=  b  sin0^.cosVgSin(4)^.-4)^) 


f  = 

Rt  r2 


fg  =  b  cos0^.cos(4i^-4^  )cosVg-a  sin0^sinvg 
sin0^sin(4)^-4)^)cosvg 

=  b  sin0^.sin(4)^-4>^)jD5gCosVg-SgSinVgj-a4)^  9^gCOS0^ 

ab4,  r  -i 

^  ^  '^Se 

g^  =  -b  ^  sin0^cos(4)^-4)^)co5Vg-a  D^gCOS©^  +  D^g  sinvg 


9r. 


ab4.rDseSinVg 


gg  =  b  cos0^sin(4»^-4)^)cosVg  +  a4)^D^^sin0^ 


sin0^cos(4)^-4)^)cos\>g 
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r?  7?  ?  ? 

=  va  sin  f  b  cos  v 
Sc  *  0  e 


?  ? 

( a  -b  )cos\)  sinv 
e  e 

r?  2  2  2 

va  sin  v +b  cos  v 
e  e 


Therefore,  one  obtains  (Av^.Aij)^)  for  a  known  ( AR^,A0^,A<t>^)  using  Equation 
(74).  In  order  to  determine  the  initial  diffraction  point  for 

a  receiver  location  (R^.,0^,(t)^) ,  one  can  always  assume  a  diffraction 
point  at  the  source  (^^,4)^)  =  (0,0)  with  the  radiation  direction  (9^,Tt/2) 
for  the  positive  ray  ( in +y  direction)or  (0^*  3n/2)  for  the  negative  ray  (in 
-y  direction),  and  gradually  add  increments  (AR^,A6^,A(t)^)  until  the 
desired  receiver  location  (R^»0^»4'^)  is  reached.  After  the  initial 
diffraction  point  is  identified  by  (Vg,i|)^);  y*  and  therefore,  the  geodesic 
path  is  determined  by  Equation  (63).  Such  a  numerical  approach  is  il¬ 
lustrated  in  Figure  0.  Note  that  one  need  not  trace  out  the  complete 
<|codcsic  path  from  the  source  location  to  the  diffraction  point  for 
each  new  radiation  direction.  As  shown  in  Figure  9,  the  diffraction 
point  (Vg+Av  g,(t)^+A(|)^)  for  the  next  receiver  location  is  determined 
from  (Vg,({i^),  using  Equation  (74),  if  (AR^,A0^,A(J»^)  are  small  which 
is  the  case  when  a  complete  radiation  pattern  is  computed. 

After  the  geodesic  path  is  determined,  various  other  parameters 
associated  with  Equations  (15)-(18)  must  be  found.  The  Fock  parameter 
of  Equation  (12)  is  given  by 


•  I  i&)' 


where  the  integral  is  evaluated  along  the  geodesic  path.  Note  that  S, 
is  the  arc  length  along  the  geodesic  and  given  by 
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s  s 

_  _  2*1 _ 

“  siriY  cosy  * 

Therefore, 


1  I  9  ?  9  9 

dg,  =  Va  sin  v+b  cos  vdv 


or. 


dil 


1 

COSY 


ad(|». 


and  the  integration  can  be  rewritten  as 


2  2 

where  pg  =  !/(< ^cos  Y ) »  and  <^,<2  are  given  in  Equations  (45) 

and  (46). 

/d'l'  (Q') 

Next,  the  ray  divergence  f actor ^  defined  as  the  change 

in  the  separation  of  adjacent  surface  rays  as  shown  in  Figure  1?.  It 
is  assumed  since  the  prolate  spheroid  will  be  long  and  slender  that 
the  ray  divergence  factor  is  unity.  Note  that  it  is  unity  for  a  cylinder. 

E.  Cone  Perturbation  for  Near  Field 

When  the  source  is  not  located  at  the  mid-section  (9^^90°,  in 
Figure  14),  the  cone  perturbation  model  is  used.  As  shown  in  Figure  131  a), 
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the  spheroid  is  modeled  by  a  perturbed  cone.  The  associated  unfolded 
surface  is  shown  in  Figure  13(b).  Ify  and  6  denote  the  angle  between 

■-  /N 

t  and  t^  at  Q*  and  Q,  respectively, it  is  seen  in  Section  III-(B)  that 
they  are  not  the  same  as  in  the  cylinder  case.  In  fact, 

B  =  ra  (79) 


where 


a  cose 

with  a'  =  acosvj^  and  tanv^  =  ^-Yfrie  ‘  "’an  i  pul  at  ion,  one 

can  show  that  ^ 

(r-Sg)cos(Y-a)  =  r  cosy  (80) 


where 


The  line  going  from  the  diffraction  point  Q  (a  cosVgCost,.,  a  cosVpSin(|)^, 
b  sinv  )  to  the  receiver  location  (R^,e^,(j)^)  coincides  with  the  geodesic 
tangent  t,  and  one  obtains 


a  sinv  cos<j)^ 

-sir»((  cos(Y-a)  -  - - - sinCy-a) 

x/a^s i n^Vg+b^cos^v  g 


R|.sine^coS(>^-a  c"SvgCOS,j)^ 
- D - 


(81) 


a  sinvgSinijij. 

cos<|)„cos(y-a)  -  ■  ■■  - ■_ - sin(y-a) 

'  !???  2 

7  a'^sin  Vg+b'^cos  Vg 


R^sine^sin(((|.-a  cos^gSin,))^ 

- D - 


(82) 
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RCE 


b  cosVgSin(Y-«)  R^cosO^-b  sinVg 

I  2  .  ?  2  2~  ^ 

Jdi  sin  Vg+b  cos  Vg 

From  Equations  (81)  and  (82),  one  obtains 

RtSin0*.sin(<|)*.-<|)-) 
cos(Y-a)  =  , 

and 

a  sinVgSin(Y-Qi)  R|.sin0^cos(4)^-(J)^)-a  cosVg 

/“? - ? - ? - 

V  a  sin  V  +b  cos 

e  e 

Equation  (85)  is  divided  by  (83)  to  give 
a  sin\)g  R^sin0^cos(<t^-4)^)-a  cosVg 

b  cosVg  R^cos0^  -  b  sinVg 

Therefore,  one  obtains 


(83) 


(84) 


(85) 


(86) 


f(Ve*Y.c*)  =  (r-Sg)cos(Y-a)-r  cosy  =  0 
9(®t»*t’V'°‘^  =  0'  cos(Y-a)-sin0^  sin((^^  -  ■^)  =  0 
h(0^,4i^,Vg,a)  =  a  sinVgCos0^+b  s1n0^cosvgcos((j»^  - 


ab 


=  0 


where 


=  ^-2|^|-  sin0|.cosvgcos(^  -  p^)+  ^  ccs0^ 


a^  2  b^  2 

+  -2  cos^Vg  +  sin^Vgr  . 

Rt  R^  J 

Using  the  increment  method  again,  one  uses 


^+1  =  ^  ^Vg^'^e  +  V®  ^  ° 


(87) 

(88) 


(89) 
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9,*1  =  g,*  g„^  AVg  +  g^Ac.  *  g^Av  *  gg^ASt  *g^^A+^  *  gi,j“t  '  “ 


'>1*1  '  •'i  *  ^*'’e  * 

In  matrix  form,  it  is  given  by 


+  hp  aR^  =  0 


'  f  f  f  1 

^»e  ™ 

AVg 

9  9  9 

\>g  «  Y 

Aa 

= 

-9i-90^AOt-g^^A4>t-9R^ARt 

h  h  0 

L  Vg  g  J 

Ay 

_-h.-hQ^A6t-h^^A4>t-hR^ARt 

(90) 


I  p  2  2  2 

/a  sin  Vg+b  cos  cos(Y-a) 


"  (^-Se)sin(Y-a) 


f  =  r  sinY  -(r-S  )sin(Y>a) 

T 


=  £Os{Y-al 


sinvgcosvg  +  ^  sinVgSinO^,cos  -  -}) 
R .  t 


^  COSVg  cose|. 


7*  IT 

=  D'sin(Y-a)+sin0^cos((|)|.  -  a)ji 


+  cos^y-g) 


-  il  cosVgSine,.sin(<|>t  -  -p 


g^  =  -D'sin(Y-a) 
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It  is  seen  that  one  can  solve  for  (AVg.Am.AY),  a  known  (AR^,A0|.,A<t>|.) , 
using  Equation  (90).  To  obtain  a  diffraction  point  (Vg,4>^)  for  a  given  re¬ 
ceiver  location  ,  one  can  always  assume  the  first  diffraction 
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point  is  at  the  source  (v,({>)  =  (v^,0)  with  the  radiation  direction 

for  positive  ray  (in  y  direction)  or  (0^=^.  negative  ray 

(in  -y  direction),  and  gradually  add  the  increments  (Ae^,A(|)^)  until 
the  final  radiation  direction  (9^,4)^)  is  reached.  Again,  one  need  not 
start  out  from  the  source  everytime,  but  obtains  the  new  solution  di¬ 
rectly  from  Equation  (90),  provided  that  the  new  receiver  location 
does  not  deviate  greatly  from  the  previous  direction. 

The  geodesic  arc  length  is  obtained  from  either  of  the  following 
equations: 


acosy  =  (r-Sg)sin3i  ,  or 
Hsiny  =  r-(r-Sg)cosa  . 


(91) 


di 


-tan(Y-a)sina  +  cosa 


and  the  Fock  paramenter  ^  is  obtained  by  integrating  along  v  or  a,  i.e., 
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The  ray  divergence  factor ,  is  again  assumed  to  be  unity. 

This  completes  the  cone  perturbation  model  used  to  compute  the  geo¬ 
desic  ray  paths  on  the  prolate  spheroid  provided  the  source  is  not  mounted 
on  the  mid-section  of  the  prolate  spheroid. 

IV.  RESULTS 

The  solutions  presented  in  the  previous  chapters  are  employed  here 
to  compute  the  near  field  radiation  patterns  for  short  monopoles  or 
slots  mounted  on  a  2XxlQX  prolate  spheroid.  To  show  the  validity  of 
the  cylinder  perturbation  solution,  the  source  is  placed  at  0^=90°  as 
shown  in  Figure  14.  To  examine  different  conical  pattern  cuts,  a  cartesian 
coordinate  system  (x',y',z')  originally  defining  the  spheroid  geometry 
in  Figure  14  is  now  rotated  into  (x,y,z)  as  shown  in  Figure  15.  Note 
that  the  new  cartesian  coordinates  are  found  by  first  rotating  about 
the  z'-axis  the  angle  (|)  and  then  about  the  y-axis  the  angle  e  .  The 
pattern  is,  then,  taken  in  the  (x,y,z)  coordinate  system  with  fixed 
and  <|)  varied. 

The  calculated  examples  are  shown  in  Figures  16-24.  The  algorithm 
boundaries  shown  in  Figure  25  are  used  in  determining  whether  one  or 
two  rays  are  used  in  the  solution.  Note  that  3 is  defined  automatically 
by  the  numerical  solution.  This  is  done  by  determining  the  caustic 
angle  in  the  elevation  pattern  ^  )  and  adding  a  few  additional  degrees 
to  that  value,  i.e.,  where  2°<A3£l0°.  One  would  expect  to 

observe  slight  discontinuities  somewhere,  because  various  numbers  of 
rays  are  included  in  different  regions. 
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Figure  16.  Near  zone  radiation  patterns  for  R^=15X  for  a  short 
monopole  mounted  at  65=90°  on  a  2XxlOX  Spheroid. 


Figure  17.  Near  zone  radiation  oatterns  for  R^=15X  for  a  short 
monopole  mounted  at  0^=60°  on  a  xlOX  spheroid. 


(d)  0j.=90°,  ♦^.=90°,  0=90°  (e)  0^=90°,  ♦^.=0°,  6=90° 

Figure  18.  Near  zone  radiation  oatterns  for  R^.=15X  for  a  short 
monopole  mounted  at  65=30°  on  a  2XxlOX  Spheroid. 


Figure  19.  Near  zone  radiation  patterns  for  R^=15X  for  an  axial  slot 
mounted  at  e  =90°  on  a  2xxlOX  spneroid. 


Figure  20.  Near  zone  radiation  patterns  for  R.=15X  for  an  axial  slot 
mounted  at  9^=60°  on  a  2XxlOA  spheroid. 


Figure  21.  Near  zone  radiation  patterns  for  R^=15X  for  an  axial  s 
mounted  at  0^=30°  on  a  3^  xlO  x  spwroid. 


Figure  22.  Near  zone  radiation  patterns  for  for  a  circumferential 

slot  mounted  at  0  =90°  on  a  2XxluX  spheroid. 


Figure  23.  Near  zone  radiation  patterns  for  R*=15X  for  a  circumferential 
slot  mounted  at  0  =60°  on  a  2Xxim  spheroid. 


gure  24.  Near  zone  radiation  patterns  for  R^.=15X  for  a  circumferential 
slot  mounted  at  6.=30‘^  on  a  2Xxnj\  spheroid. 


TANGENT  PLANE 
TO  SOURCE 


Figure  25.  Cone  boundary  used  to  define  terms  to  be  included 

in  the  shadow  region. 


V.  SUMMARY  AND  CONCLUSIONS 


The  object  of  this  study  has  been  to  develop  an  efficient  numerical 
solution  for  the  high  frequency  near  zone  radiation  patterns  of  a  spheroid- 
mounted  antenna.  The  Geometrical  Theory  of  Diffraction  is  used  in  this 
study  to  calculate  the  radiation  patterns.  The  curved  surface  diffraction 
solutions  are  discussed  in  Section  II,  where  the  creeping  wave  solutions 
in  the  shadow  region  are  of  particular  interest. 

Cylinder  and  cone  perturbation  methods  are  presented  in  Section 
III  to  simulate  the  geodesic  paths  on  a  spheroid,  which  in  turn  can 
be  used  to  model  an  aircraft  or  missile  fuselage.  Because  the  cylinder 
and  cone  are  developed  surfaces,  the  geodesic  paths  can  be  found  on 
the  unfolded  planar  surface.  For  a  given  radiation  direction  in  the 
shadow  region,  the  geodesic  path  and  final  diffraction  point  on  the 
spheroid  can,  then,  be  found  via  an  efficient  numerical  approach. 

The  Geometrical  Theory  of  Diffraction  is  combined  with  the  per¬ 
turbation  methods  in  section  IV  to  obtain  various  radiation  patterns. 

These  solutions  are  very  useful  in  predicting  the  high  frequency 
radiation  patterns  for  antennas  mounted  on  a  spheroid.  The  cylinder 
and  cone  perturbation  techniques  are  highly  efficient  in  approximating 
the  geodesic  paths  on  the  spheroid  for  this  problem.  Further,  this 
perturbation  method  can  be  extended  to  more  general  convex  surfaces 
such  as  ellipsoids. 
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